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ABSTRACT 


Title of Thesis : The Fractal Analysis of Soot Agglomerates 
Roberto J. Samson: Master of Science, 1986 


Thesis directed by: Dr. James W. Gentry 
- Professor 
~ Department of Chemical and Nuclear Engineering 


The fractal behaviour of soot agglomerates formed by the combustion 
- of acetylene in a coannular diffusion burner is studied. Structural data 
from electron micrographs were obtained by two methods: the first by 
particle counting with the aid of stereopairs for small clusters, and the 
second by electronic Hae with advanced image processing tech- 
niques, used for the larger agglomerates. An average agglomerate length 
to width ratio of about 1.7 with a standard deviation of 0.6 was ob- 
tained. The agglomerates were found to have a fractal dimension D, of 
about 1.8-1.9 based on measurements of the pair correlation function for 
large agglomerates and 1.5-1.6 based on measurements of the radius of 
gyration for small agglomerates. These results are compared with recent 


computer simulations as well as experimental studies of inorganic smokes. 
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ia Introduction 


Aerosols are suspensions of small particles in air or other gases. 
Whether formed by the conversion of gases to particulate matter or by the 
disintegration of liquids and solids, aerosols can have a strong in- 
fluence on our health and environment. One of the major challenges in 
aerosol science is ue develop an understanding of the relationships 
between structures and growth mechanisms. In regard to spherical 
particles, there has been enormous progress in the description of their 
properties and dynamics. One important aspect of growth is that due to 
coagulation, which refers to the collision of particles leading to a 
reduction of the total number concentration. For the special case of 


droplets, coagulation proceeds by the collision of two particles, which 


stick together, and finally coalesce [Friedlander,1977]. 


In the combustion of fossil fuels, ultra-fine soot particles (which 
are somewhat spherical) are produced from a complex growth process 
involving nucleation, surface growth (or condensation) and coagulation 
[Siegla,1981]. These particles may initially fuse upon collision. 
Eventually, however, the sructural integrity of the particles are 
maintained, and a cluster made up of these individual primary particles 
are formed. At the present time, our understanding of the coagulation 
process resulting in such non-spherical structures is relatively limited. 
Compared to spherical particles of the same mass, the growth rates of 
agglomerates will be more rapid due to their greater collision cross- 


section. 


Soot is generated in this manner by automotive engines, industrial 
processes, and forest fires. An electron micrograph of the soot produced 
by a laminar diffusion acetylene flame is shown in figure 1. The problem 
of size classification has often been treated by comparison to the 
size of a sphere that behaves similarly in regard to a specific property. 
One example is the use of ‘aerodynamic diameters’ in impactor measure- 
ments to obtain size distributions. Agglomerate ’diameters’ based on one 
property often differ from that based on another. It is therefore 
evident that a more fundamental description of agglomerate structure 


is needed. 


Recent simulations of particles undergoing Brownian motion with the 
constraint that the particles stick but do not coalesce upon collision 
yield clusters that are remarkably similar to actual agglomerates 
[Meakin,1984]. More importantly, the clusters were described as fractals 
[Mandelbrot,1977]. That is, a power law relationship between the mass, 
or number of particles N, and the size of the clusters R (e.g. the radius 


of gyration R,), exists in the form 


N= kk # (1) 


where D, is the fractal or Hausdorff dimension and k is a proportionality 
constant. Experimental studies of actual agglomeration processes have 
also shown similar mass-size relationships. Forrest and Witten [1979] 


obtained fractal dimensions between 1.5 and 1.7 in their analysis of Fe, 


(la) 


(1b) 


Figure 1. Samples of soot agglomerates formed by acetylene combustion. 
Masnitications: at) (a) l0Q°KX, and’ (b) 10 KX: 


Zn and SiO, smokes from transmission electron microscopy (TEM) images. 
Weitz and Huang [1984] determined gold colloid agglomerates to have a 
fractal dimension of about 1.70 from TEM and small angle neutron scatter- 
ing measurements. Feder, et al [1984] reported protein aggregates to 
have very low sticking orobapilicias resulting in more compact structures 


with a dimensionality of about 2.53. 


Medalia and Heckmann [1967,1971] characterized primary agglomerates 
(agglomerates broken down by sonic dispersion) of carbon black in terms 
of bulkiness, anisometry, and shape factors. The projected area of each 


agglomerate, A_, was measured as a function of number of particles N,- 


YS 


It was determined, for 5 < N, < 250° that 
; a 
Ne = (A, /A, ) a=1.1 (2) 


e 


where a is the projected area exponent and A, is the average area of a 
primary particle in an agglomerate. One cannot infer, however, any 
fractal behaviour on the basis of their results. In fact, for struc- 
tures with a fractal dimension less than 2.0 (as suspected for agglome- 
rates), little overlap in mass occurs upon projection, and the cluster is 
essentially transparent when viewed at any orientation. In the limit of 
large NS the fraction of particles evident upon projection becomes 


constant. As a result, N, is directly proportional to A,, and therefore a 


has a limiting value of 1.0. 


Hence the question of whether this important class of agglomerates 
exhibited fractal behavior remained to be investigated. This is the 


4 


motivation of our study, which focuses on soot agglomerates generated in 
a coannular diffusion burner. We attempt to characterize the structures 
in terms of the fractal dimension ened advanced digitization and 
processing techniques, on electron micrographs of agglomerates ranging in 
sizes as little as 50 mm to those as large as 40 um. The primary 
objective is to compare experimental and computer modelling results of 
cluster-cluster aggtomeration, as well as similar studies performed 

on inorganic smokes and gold colloids. Secondly, the effect of cluster 
shape on the accuracy of the experimental methods and future work 
involving light scattering (in terms of the pair correlation function) is 
investigated. It is hoped that the fractal information acquired may in 
the long run be applied to the study of light scattering and absorption, 


sedimentation, and diffusion of agglomerates. 


1 borer ractals 


A. Background 

In his 1977 essay ‘Fractals- Form, Chance, and Dimension’, Mandel- 
brot introduced a new geometry of nature to describe many of the irregu- 
lar shapes around us which he calls fractals. The most useful fractals 
involve chance and both their regularities and irregularities are 
statistical. In addition, the shapes described tend to be scaling, 
implying that the degree of their irregularities are identical at all 
length scales. In his analysis, the fractal or Hausdorff [Hausdorff,191- 


9) dimension forms the basis for all subsequent developments. 


The formation of clusters by the union of many small elements, 
termed agglomeration, is an extremely important phenomenon in a variety 
of processes. These include coagulation in aerosol and colloidal 
physics, percolation and nucleation in phase erane gear and critical 
phenomena, and gelation and Bile ear os processes in polymer science, 
among others. The highly branched and chain-like structure has often 
been inadequately represented as spherical in accounting for its size in 
equations of mass and momentum transfer. Witten and Sander [1983], 
suggested that aggregates of Fe, Zn, and SiO, may be classified as 
fractals. It is only recently thatwiractal concepts are being applied to 
other agglomeration processes as well. In the following review of 
fractals, we shall introduce definitions with the aid of rigid geometric 
constructs of 1 and 2 dimensional fractals, then lead to the more complex 


analysis of clusters and random points. 


B. Definition of a Fractal 
Given a set of points in space, the following relationship holds 


true’ ifsthe set isa fractal. 


_K(a) =a (3) 


where D, is the fractal dimension, and K(a) is the number of segments (in 
1 dimension) necessary to completely ‘engulf’ or cover the entire set of 
points. Similarly, K(a) for a 2 or 3 dimensional set would correspond to 
the number of circles or spheres, respectively, of radius a necessary to 


cover the set. 


In order to understand fully the interpretation of D,; by equation 3, 
consider the following fractal shapes. Although unrelated to natural 
systems, they should provide a better understanding of fractal geome- 
tries. The first set is a family of 1 dimensional curves Mandelbrot 
calls a triadic Koch curve, shown in figure 2. The points belonging 
to the set make up the intervals defined by each straight segment. 
Figures 2a, 2b, —2c and 2d ec 4, 16, and 48 intervals, respectively. 
Each succesive shape is generated by replacing the intervals’ midthird by 
the promontary of an equilateral triangle, resulting in a four-fold 
increase in line segments. The line segments have lengths one-third 
that of the previous interval. Referring to equation 3, it can be easily 


shown that the following quantities: 


figure 2 (a) K(a) 


a 1 ja 
b 1/3 4 
c 1/9 16 
d ose oe 48 


lead to a fractal dimension D,= (log 4)/(log3) ~ 1.26186. Note that K(a) 


a approaches zero, all within a confined region. 


tends to infinity as 


Next we consider a two ai heer ea Hausdorff set known as a 
Sierpinski Carpet, illustrated in figure 3. The ‘initiator’ is a 
filled-in square (figure 3a), and successive sets are generated by 
dividing each blackened square into nine equal portions and ‘removing’ 
the center square, and so on. K(a) in this set represents the number of 
squares of side a necessary to cover all of the darkened areas which 


represent all the points in the set. From the figure, we get 


figure 3 (a) K(a) 
a 1 1 
b Ws 8 
S 1/9 64 
d Wi] 514 


and D,~) 1.89: from equation 3. The Koch curve Sierpinski) carpet) and the 
other ‘designed’ fractals described by Mandelbrot all share two important 
features: the first is their self-similarity, which refers to how a part 
of the set resembles the whole, and secondly, all obey equation 3, 
which explicitly describes their scaling behaviour. Note that the 


self-similarity ratio a was not chosen at will but must be part of a 


ae 


2b 


2d 


Eirguce 2, A Triadic Koch Curves Deel 2618.7)(a) the initiator, and 
(b)-(d) successive sets generated by replacement of previous 
intervals’ midthird by the promontary of an equilateral 
triangle. Straight line distance from end points is 1.0. 


Fi curesge 


A Sierpinski Carpet, (a) initiator, 
successive sets, D,-1.89. 
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and (b)-(d) 


strict scale, namely 1, 1/8, (1/8)*, and so on. 


In contrast, fractals encountered in nature involve a certain degree 
of randomness. As a result, the self-similarity is less rigid and more 
Statistical. Self-similarity for the case of agglomerates is used in the 
sense that the relative proportions of the numbers and sizes of void 
spaces within the agglomerate are identical when viewed at different 
magnifications. Due to their finite size a limit to this self-similarity 
is approached when viewed at sizes near that of a primary particle and 
beyond the overall size of the agglomerate. The following discussion 
centers on how commonly measurable quantities such as density, number of ~ 
particles, and radius of gyration can be used in determining the fractal 


dimension. 


C. Agglomerates as Fractals 

Suppose we havea 3-D cluster made up of a set of N, individual 
particles or '’points’. Recall from equation 3 that K(a) is the number of 
spheres of radius a needed to cover the set. The average number of 


points within each sphere of radius a must then be 


Ds 
ee, - N,/K(a) ve Nia (4) 


Equation 1 is derived in the following way. Referring to equation 


4, it is obvious that as a approaches the size of the agglomerate R,, 


, — 


K(a) approaches 1, and N,,, - N,. Therefore 


ll 


Nera | (1) 
as presented. 


A new quantity called the pair correlation fuction c(r) is now 
introduced, and its significance shall be discussed shortly. By defini- 
BLovrecor) is the average density of points a distance of r from any 
point. The pair correlation function quantitatively describes the local 


density fluctuations from the overall density: 


e€r) =L1/C2xrNn) ) < p(r+e ptr) > | (33 


r 
where <...> represents orientational averaging and p(r) is the density 


at position r, and for TEM images, can be taken as 1.0 at the center of a 


particle and 0.0 elsewhere. If the set is a fractal, then 
a 
Nes | ecred’r (6) 


where d is the Euclidean dimension. It follows that an expression for 


e¢r) "in themtorm 


Se) 
m 
’ 
Q 
on 
~J 
Vw 


for a fractal set satisfies the relationship described by equation 4 when 


substituted into equation 6. Consider a random distribution of points in 
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d space which has a uniform density. Since the local density is equal to 
the overall density, D,=d. Similarly, D,=-d for solid geometric shapes 


such as spheres and circles. 


The pair correlation function has been widely used in the study of 
fluids and fluid mixtures especially in regard to the fluid properties 
near the critical point [Kirkwood and Mayer ,1951]. As the critical point 
of a liquid-vapor system is approached, molecules appear to cluster 
together in small islands. It was precisely this analogy to agglome- 
ration systems that led Witten and Sander [1981] to believe that computer 
simulation clusters of diffusion limited aggregation would have similar 
structural features. In terms of structure at the critical point, the 


pair correlation function exponent n, 


o(r) — £ 
is defined. The derivation of equation 8 did not involve fractal theory. 


Thus for systems at the critical point, the quantity 
2-7 (9) 


plays the role of D, as shown by equation 8. One must not associate 
critical point behaviour with the agglomeration process. Of interest is 
the end result of both processes, and the application of structural 


analysis techniques (in terms of c(r)) specifically to agglomerates. 


The foregoing development of the pair correlation function was 


£3 


presented in light of two main objectives in this study. The first, and 
‘most important, is it provides a means of calculating the fractal 
dimension D,. Secondly, we wish to validate these findings via light 
Scattering measurements in the near future. It was determined [Stan- 
ley, 1971} that the intensity*of light’ scattered in the k. direction, 


I(q), can be obtained by taking the fourier transform of c(r): 
1(q)/1°(q) = 1/n | d(r)e °?°*e(r) = 1/n S(q) (10) 


where q = momentum tranfer vector, defined as 


q = Ja = 2ksin(6/2) (see following figure) 


Scattering geometry used to define q [Stanley,1971] 


and I°(q) = scattering intensity for no particle correlation 
S(q) = structure factor, defined as the spatial fourier transform 
Siac un 


n = number of particles 
For isotropic particles, it was shown [Stanley,1971] that 


1(q)/I°(q) = S(q)/n~ q 727" (11) 
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provided c(r) is given by equation 8. The relationship between I(q) and 
c(r) gives an independent method for ejimet te D,. An advantage of 
light scattering is that the clusters can be studied in a medium where 
agglomeration is taking place, allowing the analysis of dynamics as well. 
In addition, the agglomerates are considered as 3-D species as opposed to 


TEM images where projections of 3-D objects are examined. 


From the pair correlation function (equation 7) , the definition of 
D- (equation 3), and the mass-size relationship (equation 1), three 
methods by which D, can be obtained is now available. In the following 
section, a review of past and present simulations of agglomeration 
systems is presented, the results of which will be compared to our 


experimental data. 
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III. Simulations 


For more than two decades, computer simulations have greatly 
contributed to our understanding of cluster growth and morphology. In 
the earliest models of agglomeration, particles were assumed to follow 
random linear trajectories. A cluster is formed by the addition of 
single particles in a random manner to a fixed particle located at the 
center of a square lattice [Vold,1963]. A later treatment included 
collisions between small clusters as well as primary particles [Suther- 
land and Goodarz-Nia,1971]. These models, however did not include the 
Brownian trajectories followed by the particles and clusters. As a 
result, the phage 34. effects which stems from the fact that a random 
Brownian trajectory has a fractal dimensionality of 2.0 (since it 
will ultimately cover an entire plane) , while a linear trajectory has a 
dimensionality of 1.0 were not reproduced by these models [Meakin,1984]}. 

Recently, Witten and Sander [1981] have developed a model that 
included the addition of single particles following a Brownian trajectory 
(through random walks) towards a growing cluster fixed at the origin. 
The resulting clusters exhibited remarkable scaling properties in 
terms of their growth and internal structure, with a fractal dimensiona- 


lity of ~ 2.5, but appearing much denser than actual agglomerates. 


For purposes of comparison to actual agglomerates, a similar 
program was created and is listed in appendix C. In this program, an 


iteration procedure is followed in which an individual particle diffuses 
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a distance of 1.0 lattice spacing from its original position, with the 
constraint that the particle stick and remain stationary upon contact 
with the growing cluster (initially a single particle at the origin). 
After contact, another particle is added at a random point along a radius 


of 3r,,,t) from(thesorigin; wherecrysypisithegdistance ¥from*the ‘origin to 


x 
the furthest particle in the cluster. The new particle again diffuses 
until it collides with the cluster. To decrease execution times, a 
particle diffusing to a distance greater than r,,,+10 is allowed to 
travel distances longer than one lattice spacing per move. In addition, 
particles¥atdistancevofv2 (ro) 945) fromethejorigintare’ considered’ lost’ 


and destroyed, and another one added at the previous starting radius 


{[Meakin,1983]. 


Although the growth from a single stationary ‘seed’ seems unrealis- 
tic, this model set a precedent for future simulations. Meakin[1984] 
proposed a more realistic model in which a fixed number of particles, 
randomly distributed inside a box (cubic lattice), are allowed to diffuse 
via random walk with the constraint that the particles stick upon 
collision. This process is termed cluster-cluster agglomeration as 
opposed to Witten and Sander’s diffusion limited aggregation or DLA model 
which limits diffusion only to individual particles. Periodic boundary 
conditions are applied to prevent loss of particles through the walls of 
the box. These simulations, carried out until one cluster remained, 


resulted in clusters with rather tenuous structures that satisfyingly 


resemble actual agglomerates, and had fractal dimensions from 1.7-1.9. 
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More recently, Mountain and Mulholland[{1985] applied concepts 
commonly used in molecular dynamics to the problem of cluster-cluster 
-agglomeration. In lieu of the previous Monte-Carlo techniques, the 
particles and clusters rigidly diffused according to the Brownian motion 


rule governed by the Langevin equation (in 1 dimension): * 
d(mvin) /dt-.=ssemivaitof (12) 


where m is the mass of the agglomerate, f the aerosol relaxation time and 
f.. a stochastic) force » The: quantity v, ‘is the velocity “inthe x dire@eion 


defined as 
dx /dt) =e v_ (439 


This allowed the history of the growth process to be studied in detail by 
straightforward comparison of the kinetics of a real system. Derived 
expressions for # allowed simulation in both the free molecular and 
continuum regimes. The solution to (14) assumes the particles are in 


thermal equilibrium with the background gas: 


r= 1, + 1/ A(v+v,)tanh(fh/2) + B, (14) 


(15) 


where r and v are the positions and velocities at time tth, given r, and 


v, for time t. B, and B, are stochastic integrals satisfying 
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<B,> = <B,> = <B,.B,> = 0 (16) 
<B,?> = 3k,T/m[1-2tanh(-28h) (Ag) 


<Soe = 6k, T/mA* [ Bh- 2tanh(gh/2) ] (18) 
where k, is Boltzmann’s constant and T is the absolute temperature. The 
reader is asked to refer to a paper [Mountain and Mulholland,1985] for a 


detailed discussion of the simulation beyond the scope of this report. 


Plots of log(R, ) versus log(k=number of particles) yielding D, 
values between 1.70 and 1.90 were obtained. Figure 4 is one such plot 
averaged over two simulations in the free molecular regime for fr = 0.05 
( r is the time required to free stream a particle diameter, set = 1.0) 
and density = .0167. Greater scatter for large R, is observed due to an 
inherently smaller number of the large relative to the small agglomerates 
produced ™s These results (in terms of De) agree quite well with ‘Meakin“s 


Monte-Carlo simulation. 


a) 
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Figure 4. Ro vs- an plot for Mulholland simulation in the free molecular 
regime, five run average. ; 
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IV. Experimental Procedure 


A. Generation of Soot 
1. Overview 
Figure 5 shows the apparatus used to produced the soot agglomerates. 
The three principal components of the system are: 
(1) a coannular diffusion burner, 
(2) a glass chimney and dilution chamber, and 


(3) flow controllers. 


The burner design is similar to that used by Santoro([1983] to obtain 
particle concentration profiles within the flame. Acetylene (99.9%), is 
chosen as the fuel source due to its high sooting tendencies. The fuel 
flows through a 1.27 cm (1/2 inch) nominal brass tube at the center of 
the burner. Air is introduced through two inlet ports at the bottom at a 
rate of 40 1/min and flows laminarly and uniformly across the top. 
Pre-mixing of air and fuel is not allowed prior to combustion (diffusion 
flame). Once the proper air and fuel flow rates have been selected, the 
fuel is ignited, and a glass chimney is placed over the burner. A fine 
thread of smoke emerges from the laminar flame which extends to a 

height 30-35 cm from the base of the flame. At this height some back- 
mixing occurs due to the flow constriction at the top, as evidenced by 
the dispersion of the fine thread of smoke in this area. Agglomeration 
begins in the flame and proceeds as the clusters traverse the length of 


the thread. Large agglomerates, in the order of about .2-.5 mm, precipi- 
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Figure 5. Schematic of soot generation apparatus. 
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tate out of the region below the funneled top at high fuel flow rates, 
suggesting that agglomeration occurs to some degree by backmixing. The 
species are sampled by briefly passing TEM carbon coated grids” (3 mm 
diameter, .13 mm thick, and 200 mesh copper) over the top of the column. 
We believe that the agglomerates are drawn towards the grids by themopho- 
resis [Davies,1966] as opposed to impaction or diffusion. The relatively 
weak dependence of Parner cir on particle size allows size ranges 
between three orders of magnitude to be collected [Hinds,1982]. Samples 
were collected only for a short period of time to avoid impaction of 
clusters on top of one another on the TEM grid. A previous Pe tied ce 
and Huang,1984] of gold colloid agglomerates produced in aqueous solution 
indicated a flattening and hence distortion of the clusters on the TEM 
grid during the sample drying/preparation process. In our investigation, 
these procedures were not required prior to TEM analysis. Thus the 3-D 
structural features are preserved, as later SOS examination of the 


clusters from different viewing angles. 


2. Burner Design and Operation 

A schematic of the coannular diffusion burner is given in figure 6. 
The design specifications require that the flow of acetylene should be 
laminar and variable between 40 and 150 cm?/min. To insure this a 1/2" 
nominal brass tube (pure copper is reactive to acetylene) is used, 
limiting the Reynolds number to < 10.0 at the maximum flow of 150 
cm?/min. A tube length of 16 cm eliminates possible entrance effects 


and assures fully developed flow (length > .04 Re). 
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In order to maintain a stable flame, assure a sufficient supply of 
oxygen for combustion, and to confine the smoke produced into a finer 
thread above the flame, a laminar stream of air flows vertically through 
the top of the burner. As noted earlier, air is introduced at the 
bottom through two inlet ports. The air is dispersed uniformly inside 
the burner as it flows through three chambers, essentially composed of 2 
empty chambers betwéen two meshed screens a distance of 3.8 cm apart, 
packed with 1/4 in.(.64 cm) diameter glass beads. Laminar flow is 
achieved by passage through a 1 in. thich ceramic honeycomb with 1 mm 
vertical channels. As also mentioned, a glass chimney with a height of 
53 cm (about .3 cm thick) encloses the entire burner from its base to 
about 40 cm above it. This proved to be crucial in preventing external 
air currents caused by laboratory exhausts from distorting the steady 
flame and soot path. In addition, the funnel shaped top provides an area 
of greater soot concentration suitable for sampling. The stack was 


constructed of glass to allow the operator to visually monitor flame and 


sooting conditions. 


B. Soot Characterization 

Cascade impactor and sooting efficiency measurements are taken by 
sampling a fraction of a cross-sectional path followed by the aerosol. 
The soot concentration must therefore be uniform across the sampling 
area. To meet this criteria, the gas and particle mixture exiting the 
chimney is diluted with an additional 30 liters/min of air and is. 


dispersed through a 1 inch hole centered on a rectangular aluminum 


tripper plate. A second stack (mixing chamber) placed directly on the 
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Figure 6. Burner schematic showing glass beads location and ceramic 
honeycomb. Plate at base of burner supports glass chimney. 
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plate provides a path where turbulent flow is allowed to develop, and a 
.64 cm (1/4 inch) glass tube situated at the center of the stack 30.8 cm 
above the plate is used for sampling. We note that once the soot is 
dispersed, futher agglomeration is unlikely because of increased inter- 
particle distances. 

Size distributions for two flow rates were determined using an 
Anderson 2000 cascade impactor. The impactor is comprised of eight 
stages, each Beet vere of an aluminum orificice plate and a base for 
filter support (figure 7). The stages are sealed with O-rings and held | 
together by three spring clamps. Orifice diameters are the same on each 
plate, but get increasingly smaller in each succeeding stage. When air 
is drawn through the sample, impaction on any given stage depends on the 
particle aerodynamic diameter, defined as the diameter of a hypothetical 
sphere of unit density with the same Stokes number (or settling 
velocity) as the particle in question [Friedlander,1977]. The jet 
velocity of the stage and the cut-off of the previous stage determines 
the range of particle sizes collected on each stage. Particles not 
collected on one stage follows the gas stream around the edge of the 
plate to the next stage, where it is impacted or passed to the next, 
until the jet velocity is sufficient for impaction. The air from the 
burner is sampled at 1.0 cfm for a pre-determined period, depending on 
particulate concentration. At this flow the particulate fractionation 
range is limited from .4 to 10 wm diameter. Figure 8a and 8b are size 
distributions at fuel flows of 42 and 69 cm?/min, respectively. The 


average diameter is shown to increase by a factor of ten between the two 
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Figure 7. Andersen 2000 cascade impactor diagram [Andersen,1977]. 
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Figure 8a. Size distributions based on aerodynamic diameter from 
impactor measurements at fuel flow of 42 cm? /min. 
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Figure 8b. Size distributions based on aerodynamic diameter from 


3 ° 
impactor measurements at fuel flow of 69 cm’ /min. 
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flow rates, from .48 to 5.8 um. 


The flows of air and fuel must be steady to insure consistent 
measurements. Air flows are controlled by Tylan FC 261 and 262 mass flow 
meters that adequately compensate for line pressure changes. A Matheson 
602 ball meter (calibrated with a bubble meter) was used to vary fuel 
flow. Table 1 is a list of sooting efficiencies at different flow 


races: 


Acetylene Flow (cm?/s) mg soot produced per 
mg C,H,_burned 
0..92 .0445 
TwLo .0639 
1.46 .0671 
Pia Vas! .0699 


Table 1. Soot production at various flow rates. 


Samples ranging from 0.1 to 1 wm were determined to be best obtained at a 
fuel flow rate of ~ 1.0 cm’/s . Larger agglomerates (5-38 yum) were 
sampled at a rate of ~ 1.5 cm’/s. Prior to TEM analysis, preliminary 
specimens where collected on microscope cover slips for viewing on a 90X 
objective lens to determine appropriate flow conditions. TEM micro- 


graph digitization procedures will now be discussed. 


C. TEM Image Processing 
The TEM images can be classified into two cluster size categories. 
The first category make up clusters about 200 primary particles or less, 
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to be used in the ‘'particle counting’ method. Magnifications in this 
category ranged from 73 to 330 KX, depending upon cluster size and image 
resolution. In this method, the micrographs were enlarged to 8 x 10 
photographs, and the center of each primary particle recorded on an 
Hewlett Packard 87 with a digitizing tablet. Figure 9 is one such image. 
Toyaia vin distinguishing overlapping particles, stereoscopic pairs of the 
images were also used. To create a stereoscopic yake. two electron 
micrographs of the same cluster are made, with the second image taken 
after the sample is rotated approximately 10° from the axis of tilt. The 
Jeol 200cx transmission electron microscope has a viewing stage that is 
capable of tilting the TEM grid (and thus the sample) by 30° in either 
direction. When viewed side by side through special lenses, one is able 
to perceive three dimensional features of the cluster. An example of a 
stereopair is shown in figure 10. Gray, et al[{1985] found that the 
number of particles counted in large samples of diesel particles was 
about 36% greater using stereoscopic pairs as opposed to conventional 


viewing of 2-D images. 


Agglomerates between 5 and 40 um constitute our second size catego- 
ry. In this set, each micrograph is converted to an electronic grid of 
512 x 512 pixels. Each pixel, as viewed on a high resolution monitor, 
has an associated gray level which can have values from 0 to 255. A gray 
level of O signifies the lightest gray, while 255 the darkest gray level. 
Using an image processing program known as Lispix, appropriate gray level 
‘thresholds’ are chosen to determine which pixels represent areas on the 


grid that are occupied by a particle. For example, if a digitized image 
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AGGLOMERATE SIZE PARAMETERS 


2RAg/d=14.2 
(LW) 2 /d=19.2 


Figure 9. Digitized micrograph of a cluster consisting of 164 particles, 
illustrating size parameters L, W, and R,. 
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Stereopair of the agglomerate shown in figure 9. 


Figure 10. 


a3 


is assigned a value of 155, all pixels with levels from 155 to 255 are 
reset to 255 (white), while pixel gray levels from 0 and 154 are set to 0 
(black). We note that the micrographs are negatives whereby clusters are 
light gray set against a dark background. The end result is a 512 x 512 
array representing the threshold version of the digitized image. For 
simplicity, values of 255 are set to 1.0 so the array consists of only 
I’s and Odds. Aa avee ites of 1.0 are termed ‘occupied sites’ since thay 
signify the presence of a particle, in either a portion or its entirety, 
at that location. Surrounding clusters not associated with the cluster 
of interest are eliminated from the threshold version by assigning 
corresponding pixels a value of 0.0 through Lispix. Data in this 


category are evaluated by the successive square, pair correlation, and 


Sullivan [1985] methods. 
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V. Data Analysis 


In this“section, the techniques*used to determine D,; will be 
discussed in greater detail. We noted in the previous section that all 
the data pertaining to the D, were obtained from TEM micrographs of 
agglomerates. Consequently, it 1S sour Pee to obtain 3-D structural 
information from the 2-D TEM images. The data analysis procedures to be 
discussed later in this section depends on the assumption that essential- 
ly all particles are evident in the projection, i.e. that the agglomerate 
is transparent and only a few particles are completely masked by others. 
The validity of this assumption relies on the fractal dimensionality 
being less than 2.0, which is believed to be the case for agglomerates. 
Mandelbrot{1977] states that D, < 2.0 implies that the projection 
typically itself is of dimension D,. Weitz and Huang[1984] offers a 


simple argument based on the relation 


MEet Re <2 (19 


which is the same form as equation 1 but N is replaced by the cluster 
mass M. Upon projection, M does not change, while in two dimensions 
normal to the projection, R also remains constant. Thus the projection 
remains a fractal, with the same D,, provided the mass can be measured. 
This will be possible only if there is not too much overlap in the mass, 
which produces the condition D, < 2. We shall further discuss the 
validity of this statement in light of the analysis of simulation 


clusters in 3-D and their projections, in addition to the previous 


a5 


verification of TEM results by small angle neutron scattering measure- 


ments in a study of gold colloid agglomeration. 


A. Particle Counting Method 

In this method, the total number of primary particles are deter- 
mined as a function of cluster Terk The clusters chosen are limited to 
oes consisting of 200 particles or less, beyond which one encounters 
increasing difficulty in distinguishing overlapping and partially fused 
spheres, even with stereoscopic pairs. Two definitions of cluster 
size are used (see figure 9): the first is in terms of the geometric mean 
of the cluster length L, and W, the projected length perpendicular to the 


axis through L, and the second definition in terms of the radius of 


gyration R,, as given by 


N 
mR, * a ) (m,xr,*) (20) 


i=l 


where m, is the total mass of the cluster and m, is the mass of particle 
i whose center. is at a distance r, from the centroid. It was assumed 
that all m,'s were equal, making m, directly proportional to the total 
number of particles N. However, particle diameters actually varied by as 
much as 16% within a cluster. In an attempt to minimize the error 
introduced by this assumption and also to account for different magnifi- 
cations, R,, as measured from the photograph, was normalized by dividing 
its value by the average particle diameter d. Values of d were obtained 


by arbitrarily choosing 10 representative primary particles from each 
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cluster. The geometric mean (LW)!/* was normalized in the same fashion 
as R,. A useful quantity is the ratio L/W, which provides a description 
of overall cluster shape. To test whether the cluster growth is fractal, 


we plotted log(size) versus log(N). For fractals a linear dependence 


with a slope of Dig will be observed (see equation 1). 
B. Analysis of Large Clusters (5-40 um) 


1. Successive Squares Method 

As previously mentioned the digitized data from electron micro- 
graphs of large clusters consisted of 512 x 512 arrays with values of l's 
(occupied site) or O's (unoccupied site). The first method to evaluate 
D, from this data set ipo ead placing series of non-overlapping squares 
of side 2£ from the center of mass of the agglomerate. The total number 
of sites occupied within each square is then determined as a function of 
£. This is similar tothe Nested Square method of Forrest-and Witten 
[1979], except that the center of all the squares are the same in 
our analysis. In the Nested Squares method, the square center location 
was chosen as the point where the center of mass of occupied sites within 
the square and the geometric center of the square coincided. We found no 
theoretical justification for doing the same. For a fractal, N will 
vary as ia A weighted fit of log(N) versus log(£) was perfomed, 
assuming that the square of the probable error for each point varied as 
N'1/2. A listing of the successive squares computer code is given in 


Appendix A. 


ay, 


2. The Pair Correlation Method 

In this method, the pair correlation approach was used to find D,. 
By definition, c(r) is the average density of occupied points in the 
digital image at a distance of r from each occupied point: The use of a 
512 x 512 lattice, however, places a limit on the number of points 
(occupied or unoccupied) that are a distance of r around each point. 
Consequently, we approximate this density as the density of occupied 
sites in a shell of inner radius r-0.5 and outer radius r+0.5. From 
equation 5, c(r) is determined by averaging the density at r over all 
occupied sites. This is done for a reasonable range of r values (1-200 
pixel spacings). For a fractal, c(r) should vary as ad , where d=2. 
since the clusters are of finite size, it is also expected that equation 
5 will be valid for a short range of r, since for longer r values, larger 


areas outside the edge of the agglomerate are encountered. 


A computer program is used to compute c(r) (Appendix B). In the 
program, the distances between each and every particle, rounded-off to 
the nearest integer, is determined. Then the number of occurences of 
each distance is recorded in array n(i). If n(51)=2100, for example, the 
total number of times when an occupied site is at a distance of r ~ 5l to 
another occupied site is 2100. The value of c(r) is computed by dividing 
n(i) by the total number of occupied sites N sae os the total number of 
points (occupied or unoccupied) in a shell of inner radius r-0.5 and 


outer radius r+.5 in a square lattice. 
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3. The Sullivan Approach 


The third and final method used for the 512 x 512 data set is a 
program developed by Sullivan({1985]. It makes use of a tree method 
which efficiently stores important structural information. This procedu- 
re can be best described by considering a 1 dimensional fractal set (i.e. 


N points scattered along a straight path) between points X, and x,: 


The set is normalized by scaling the minimum (x,) to zero and the maximum 
sUe LOLs we nevelore,, all the points belonging to the. set are represe- 


nted by a fraction between O and 1: 


The point locations are then converted to binary fractions (i.e. 1's and 
O's). For example, if x=.25, the corresponding binary fraction is ‘Ol’ 
cancers o=.0*(1/2)>+71%*(1/2)* 5 andex=./3 has a binary fraction of ‘11’ 
Sauce ./5>= 1*(1/2)* + 1*(1/2)*. These binary fractions represent 
directions or paths each point follows when a hypothetical tree is 


constructed. 


Consider the tree depicted in the following figure: 
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level 0 a=(1/2)° 
level”.1 -a={1/2)> 
level 2 a=(1/2)? 
level 3” a=ti/7) 
Starting at level 0, each point has a unique path to follow which is 
defined by its binary fraction. The binary sequence is read from left to 
right, designating a value of 1 as a direction downward to the right, and 
O as a direction downward to the left. Taking x = ‘01’, for example, the 
path 
level 0 


ee leven as 


LJ [sig Ei level 2 


is followed (shown by the dark arrows). The procedure is repeated for 
each point, while storing the number of ‘hits’ to the boxes in each 


level. At level 2, there are 4 boxes: 


[im [] [ae (| level 2 
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The number of hits contain information on how the set is structured. The 
hits in box a, for example, represent the number of points between 0 and 


0.25, in boxib, between...25°and 0.5, in boxic, between 0.5 and 0.75, and 


in box d, between 0.75 and 1.0. 


Recall equation 3, 


(3) 


where K(a) is the number of segments (in 1 dimension) of length a needed 
to cover wate set. The tree, as constructed, actually represents each 
level as a= Gavin "REL i The number of boxes that have been hit at 
each level is K(a). Thus the fractal dimension can be computed from the 
slope of log(K(a)) versus log(a). The program computes D, from K(a) 


end a atveach level. The error in D; decreases as the number of boxes in 


each level approaches the total number of occupied sites N. 


The Eocene ep provides an efficient algorithm that has a short 
execution time since K(a) is computed at values of uniform distribution 
in a log-log plot, and is suitable for vector programming on the Control 
Data Corporation Cyber 205 supercomputer. The actual program used is a 
similar version but is applied to two dimensional fractals. The x and y 
coordinates are scaled from 0 and 1, but the tree is instead constructed 


of 4 possible branches per box, i.e. for the x and y binary: 
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level 0O 
level 1 


The pair correlation and Sullivan’s method were verified using a two 
dimensional fractal known as a Sierpinski carpet with a fractal dimension 
De ~ Deo aeee eeeeion II), along with compact ordinary objects with 
D,=d=2. The fractal dimensions computed were within 2% of the true 
values,.»In) the) next, section,: the, experimentalevalues; of sD-wfor) the soot 
agglomerates, along with other structural features are presented. These 
results are then compared with simulations and past experimental studies 


of inorganic smokes and gold colloids. 
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VL. Results 


A. Particle Counting Method 

Thirty seven clusters were collected and evaluated by the particle 
counting method. Table 2 is a listing of cluster size as a function of N 
for a few of the samples collected. Our own ability to distinguish 
individual particles limited the size range to 5 < N < 200, even with 
stereopairs. The results are summarized in figure 1l. The linearity 
observed indicates a mass-size power law relationship described by 
equation’? - Fractal dimensions of 1 lol t. 913 .and ‘1.49 + .16,for cluster 
sizes defined by (LW)*’* and R,, respectively, were obtained by linear 
regression. A rather diverse set of shapes (L/W), from Hoes of high 
sphericity (L/W-1.0) to those approaching linearity (L/W>2) were obser- 


ved. The average agglomerate length to width ratio was 1.68 (a = 0.50). 


In view of our expectation that the number of particles counted may 
be low for large agglomerates due to the obscurity of particles, it is 
reasonable to expect the slope of both experimental lines to be slightly 
lower than theoretically expected for a fractal. As seen in figure 4, 
there are certain regions of high overlap where undercounting is likely. 
Below we show how undercounting of N will further reduce the value of D, 


when using R, as cluster size. 


We let n represent the difference between the actual number of 


spheres, N, and the measured number, N,. Expressing equation (20) in 


terms of N, N_, and n, we obtain 


43 


IZ 


nm 


MPOoOrDOrHrHDEFWAWONNDKH UMW WOW 


_2 (nm) Re/d (LW)L2 7/4 L/w> 


ges) a, Pa Be, ee 
4.1 2* 36 Se Gan" P 
fae Ui 4.46 eed 
Sm led 8 4.24 1 OD 
220 1.8) iad Be 
2.0 1.41 4.97 1.24 
20 Hee pie ys Yt 
bie 205 6253 ds 02 
be Eas 65.33 1.88 
2e7 2 vk 8.68 1.58 
a eee vou 1.08 
2.0 218.8 a 1.08 
2.4 bees 9.06 bea ¥ 
34,0 bom? 2 12s : Fe 
EAL 4.82 12.4 1.74 
27 4.62 La 1.68 
8.8 be ya) 14.0 ieU 
eRe) 6.24 14.3 Zr 
oes ghd 22.4 Loc e 
3.9 J else LO sez 2.06 


Typical results for the particle counting 
method. 


o refers to the standard deviation of mean d. 


average L/W for entire set is 1.68, o =.50 
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Particle Counting Method 


Experimental Data 


10 


1 


N [number of particles] 


10° 


10 10 


cluster size 


Figure 1l. R(size) vs N for experimental data. Cluster sizes defined as 


R, or (LW)?/2, 
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n 
ONRG ?oeg NOR eee ) re (21) 


where R,,, is the measured radius of gyration and r, now represents the 
distance from the center of mass to the it® sphere (i=l to n) not 


accounted for in the measurement. Rearranging equation (21), 


n Lee 
Re /R ee ae ee ) Sead gee (22) 


A two-term Taylor series expansion is performed to obtain a simplified 


expression adequate for making qualitative comparisons: 


nN 
eS mms Oe (ar ) r,7)/Ryp?-n (23) 


In general, we observed that overlaps are present near the agglomerée ze 
center of mass. Now consider a situation where 20% undercounting hes 
occured, i.e. N= 0:8N° and n= 0.2N. For the extreme case of rosv 


(overlaps at the center of mass) we get, from equation 21, 


Ree (24) 


gm 
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and ior ry ~ O.5R, , 


KR omael20sh. CED) 
Equations 22 and 23 show that undercounting further aggravates the slope 
D, by yielding superficially higher R,, values. It is therefore expec- 
ted, ~as7supportediby figure 12, 0that a lower slope of°D, will resuit 
using R, as cluster size as opposed to (LW) ET? The errors introduced 
will be more pronounced for elongated structures, typical of those 

which where evaluated. On the other hand, note that the apactred (LWyt hs 


is independent of the fraction of particles undercounted. 


B. Pair Correlation, Sullivan, and Successive Squares Method 

In table 3 , the results of the pair correlation, Sullivan, and 
successive squares metMods are listed for the size range 5.5 to 38 um. 
Agglomerates quite similar in shape to those of the particle counting 
method were collected (i.e. elongated and spherical) but on a much ele 
scale. The average length to width ratio of these agglomerates was 1.8 
(o=0.66). Table 3 shows that from each of the three methods, the 


agglomerates have fractal dimensions between 1.75 and 1.95. 


The pair correlation function showed the greatest sensitivity 
towards the internal structure and overall shape of the agglomerates. 
The power law dependence of c(r) on r clearly suggests fractal behaviour 


in terms of the scale invariant internal structure of the agglomerate. 
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However, that range of r where equation (7) is valid appears to be 
strongly dependent on the shape (L/W). Figures 12a and 12b are examples 
of clusters with high sphericity (L/w = 1.1 and 1.2) where power law is 
observed over a long range of r (2-37 and 2-24). On the other hand, 
cluster number 2513 (L/W = 2.9) shown in figure 12c is at best linear for 
the range r=2-7. This is probably due to the fact that for elongated 
structures a greater percentage of the area around the edges are encoun- 
tered as opposed to rounded ones. These areas of zero density signifi- 
cantly lower c(r) as seen in all cases as r approaches the size of the 
agglomerate. By the same token, eater agglomerates (< 5 uwm,) such as 
cluster number 2500 shown in figure 12d whose internal structures are not 
yet fully developed also display shorter power law ranges. Weitz and 
Huang[1985] compensated for this effect by considering occupied sites 
only a short distance from the center of mass. A summary of the effects 
of cluster size and shape is illustrated in figure 13a, a combined plot 
of c(r) versus r for rel to 250. Those agglomerates with L/W~l are shown 
to have longer power law dependence which is followed by a steeper drop 
for large r when edge effects are encountered. Elongated or underdevea 
loped structures followed a gradual decline towards the limiting value of 
c(r)=0. As a test case, c(r) was determined for a solid square (L/w=1.0) 
and rectangle (L/W=3.0) of similar sizes. The result is plotted in 
figure 13b. For slopes equal to 1.98 + .01, a longer range of linearity 
is found for the square. The long range effects of shape on c(r) may not 
have any direct bearing on D,. This information, however, may prove to 
be useful in any future work regarding light scattering of particles of 


such shapes and sizes. 
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table 3. continued 
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Figure 13. Long range effects of shape on c(r). (a) combined c(r) plots 
for clusters in previous figures, and (b) test case for 
filled Square (L/W=1) fit r in range 1-18, Df=1.98 + 0.01. 
and filled rectangle (L/W=3.0), D,=1.98 + 0.01, fit r in 
range 1-15. 
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There is very good agreement between the pair correlation and 
Sullivan approach for clusters between 5 and 12 wm, corresponding to 
magnifications between 5 ant 10 KX, respectively. Below 5 KX, however, 
D, by the c(r) approach becomes consistently higher (by ~.05) than those 
of greater magnification. It appears that the digitized images become 
denser for agglomerates larger than about 12 wm, most of which are 
photograped at magnifications less than 5 KX. First we consider the 
possibility jthat the increased D, for agglomerates of sizes > ~12 ym be 
the result of instrument limitations or events unrelated to the growth 
phenomena. More specifically, there is an inherent decrease in resolu- 
tion on a 512 x 512 lattice as more than one particle begin to occupy 
each pixel. Table 4 shows the approximate length between pixels in 
comparison with the average particle Sue at different magnifica- 
tions. At about a magnification of 5 KX, the length of one side of an 
occupied site is about equal to a particle diameter. The resolution 
further diminishes at 2 KX where void spaces between the particles must 
have areas greater than that projected by ~2.4 neighboring particles to 
be distinguished as occupied sites. For 1.5 KX , this area increases to 


3.2 neighboring particles. 


To more accurately measure the effects of decreasing resolution, 
the 512 x 512 array representing cluster number 2313 was reduced to a 256 
x 256 array of 1’s and 0’s. This was done by checking groups of 2 x 2 
pixels in the 51D 512 array for = 50% occupancy (values of 1’s). If 2 
50% occupied, the corresponding location in the 256 x 256 matrix is set 


to 1.0, otherwise it is 0. In effect, a 256 x 256imatvixers -obtained 
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TEM Magnification number of pixel lengths equalling 
one particle diameter (30nm) 
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Table 4. -Resolution of mass in digitization process. 
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just as if the micrograph was digitized with such a resolution. The 
procedure is repeated to obtain a 128 x 128 and 64 x 64 matrix. The data 
was then evaluated the c(r) method. As shown in figure 14, the power law 
slope is essentially the same for all four curves. Although deviations 
from power law behaviour appear at different r values, these would 
actually correspond to approximately the same lengths if plotted versus 
micrometers. We therefore conclude that the effect of low resolution 
does not result in a ‘superficially higher D> *as"suggested. ~ Ii rTetioe 
spect, if the agglomerates are indeed fractals, the internal structure is 
statistically scale invariant and changes in length scales (varying 
resolution) would unaffect D,..' Whether bigger clusters result in larger 
exponents, which implies that some other physical mechanism takes over, 
will be discussed in the next section in light of computer simulations of 


cluster-cluster agglomeration. 


To examine the poSsible effects of the imaging process on D,, 
variations in threshold levels and tilt angles were investigated. 
Micrograph number 2313 was processed at threshold levels 120, 125, and 
130, and only a 1% change in D, was found by the pair correlation and 
Sullivan methods. In general, a threshold value within + 5 gray levels 
of those that represent a contiguous image adequatély resembling a 
silhoutte of the agglomerate is chosen. In terms of the tilt angle, no 
significant effect was also found. Cluster number 2313 was tilted 55° 
and is represented by cluster 2315. As an aside, one can infer that the 
agglomerate extends out in all three directions as opposed to two. If 


growth proceeded along a plane, one would expect a prominent increase or 
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decrease in size as the sample is rotated, which is certainly not the 
case for clusters 2313 and 2516 upon rotation. This also upholds the 
assumption that the sample is not distorted or flattened on the TEM Beet 
which was the case for gold colloids when dried during sample prepara- 
tion. These results indicate a moderate insensitivity towards thresholds 
and tilt angles. Consequently, the analyst has Greste tation: ste and 


control in the digitization process. 


Generally higher values of D, resulted from the successive squares 
method. Figure 15 is a plot of log (N) versus log(2) with a slope of - 
D-=1.83. This method, however, was found to be ineffective for eluseeee 
where the center of mass is in a region of unoccupied sites, characte- 
ristic of small agglomerates with few branches, resulting in non-lineari- 
ty. In regard to computer execution times, Sullivan’s approach was found 
to be at least two orders of magnitude faster than the c(r) method. The 
advantage of the c(r) method, however, becomes apparent when insight into 


the contribution of individual particles to total scattering is desired. 
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Figure 15. Successive Squares method plot, Fit for #=2-181, slope 
D,=1.82. 
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Vile sDpisciussi0n 

The results given in the previous section now permits a qualitative 
and quantitative look at actual and predicted properties of agglomerate 
growth and morphology. Figure 16 shows comparative sizes of clusters 
obtained experimentally and those generated by DLA and cluster-cluster 
agglomeration. Note the similarities between the open and chain-like 
structures of experimental and cluster-cluster agglomerates. DLA 
aggregates appear more dense and compact, characteristic of its high 
dimensionality of 2.50. On the basis of appearance alone, it is obvious 


that cluster-cluster agglomeration best represents the soot agglomeration 


process. 


For a quantitative comparison, we applied the D, methods of 
analysis to the simulation clusters. We first address the question of 
whether D, is indeed invariant upon projection by applying the particle 
counting and pair correlation method to clusters generated from Mountain 
and Mulholland’s simulation. The availability of their program provided 
the opportunity to vary operating parameters such as the aerosol relaxa- 
tion time and particle density. Table 5 shows the results of the 
analysis for Br=.05 and p=.0167. For the particle counting method, the 
results indicate a 3-4% decrease in D, when evaluating projections 
instead of actual 3-D coordinates. Accounting for the errors in the 
slope of the N vs. R, plot, however, makes this difference relatively 
insignificant. To test the pair correlation approach, a ‘digitized’ 


projected image was created from a Monte-Carlo integration technique 
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Table 5. Comparison of 3-D and projected values of D, from cluster- 
cluster simulation. 
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listed in Appendix D. As seen in table 5, there is good agreement 


between D- in two and three dimensions. 


Experimentally, we obtained a D, of 1.49 + .13 and 1.61 + .08 by 
the particle counting method for the set of agglomerates in the range 5 < 
N < 200 for cluster sizes defined as R, and (LW)?/2, respectively. 
Large experimental clusters from 5 to 15 um were stated to have a 
dimensionality D, ~ 1.85 by the pair correlation and Sullivan methods. 
No corresponding increase in D, was experienced for the diffusion limited 
aggregation model. In DLA, as, example, the pair correlation analysis of 
a cluster of ~3000 particles yielded similar De values (~2.4) as that 
obtained for an R, vs k plot for 1 < k < 200 particles. Im contrast, 
cluster-cluster agglomerates exhibited a dimensionality of 1, 70° fen ae 
N s 200 and 2.00 for 2 < N < 2000 by the particle counting method as 
previously mentioned. This apparently suggests an increase in D, as 
observed experimentally. “However, there is greater scatter in the data 
for large N since fewer of these clusters are formed. This result is 
still therefore inconclusive. Figure 17 is a comparison of experimental 
and simulation data for the particle counting method. The much steeper 
slope for DLA is obviously the result of its denser structure. The 
fitted equations for both simulations are also valid for k up to about 


200 particles. 


It was also reported that clusters > 12 um had slightly larger D, 
values (~1.9). A decrease in digitization resolution was dismissed as an 


explanation for this finding. It is therefore probable that either 
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LAMINAR DIFFUSION FLAME WITH RESULTS 
OF COMPUTER SIMULATION 
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Figure 16. Qualitative comparison of clusters produced experimentally 
and by computer simulations. 
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Figure 1/7. R, vs N plots of experimental and simulation data. 


66 


changes in growth mechanism occur for the extremely large clusters, or 
that they could not withstand external forces put upon them in the sample 
generation and preparation process. Assuming the growth process proceeds 
by some form of cluster-cluster agglomeration, single contact points 
which link clusters may not be sufficiently strong to maintain the 
agglomerate’s overall structure. As a consequence, reshuffling of 
clusters to form “verre more stable agglomerates become likely during 
harsh collisions in backmixing or impaction on the TEM grid. To a lesser 
extent, dense structures can also arise from impaction of small clusters 
onto the sample immediately after it is collected, particularly those 
with a large surface area. Note the large number of clusters surrounding 
cluster number 2506 in figure 12b. Due to execution time restrictions, 
starting concentrations > 8,000 particles has yet to be performed for 
Mountain and Mulholland’s simulation. Meakin[1984], reports D, ~ 1.80 
for clusters made up of 30,000 particles with the Monte-Carlo cluster- 


cluster agglomeration model. 


We now compare our results to experimental values reported in the 
literature. Forrest and Witten[1979] produced Fe and Zn agglomerates by 
vaporization from an eletroplated tungsten filament. From pair correla- 
tion analyses of TEM images, the agglomerates were determined to have a 
Dy-lso-1./7.) It was reported, however, that the agglomerates were not 
examined in their entirety. The extreme sizes allowed only a fraction of 
the total area of the agglomerate to be photographed. In any event, 
their pioneering work provided preliminary evidence of the fractal 


nature of these agglomerates. 
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Weitz and Huang[1985] presented an in-depth study of the kinetics 
and morphology of gold colloid agglomerates in aqueous solutions. 
Negatively charged gold spheres (d,~14.5nm) which resist agglomeration 
due to repulsive forces were produced by the reduction of gold salt 
Na(AuCl,) with sodium nitrate. The addition of pyridine neutralizes 
the particles which then allow agglomeration to proceed. It was noted 
that the process is highly irreversible, owing to their ability to 
withstand physical or chemical separation measures. The agglomerates 
were found to have a dimensionality of 1.70 + 0.1 from structural 
' analysis of TEM projections by k versus (LW)?/2 and c(r) versus r plots. 
As with the metal smokes of Forrest and Witten, no apparent change in D, 
for small terene were reported. The study however, satisfyingly showed 
that D, is indeed invariant upon projection through confirmation of their 
results with small angle neutron scattering measurements in three 


~ 


dimensions. 


In summary, experimental values of D, ~1.80 for soot agglomerates 
ranging from 5-15 um were found to be in agreement with Brownian motion 
simulations of cluster-cluster agglomeration using Monte-Carlo and 
Langevin trajectory approaches. There are obvious differences -in results 
for small clusters (2-200 particles) and extremely large clusters ( > 12 
um). This should not imply that the methods of analysis or the simula- 
tions are totally inadequate. We believe that some physical effects are 
responsible for the transition from D,~1.5 to 1.9. Simulations currently 


take into account only a few of the controlling factors in a real system, 
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and future work in this area should emphasize possible rearrangements due 
to weak bonding resulting from single contact points and more accurate 
simulations of the contributions of rotational motion. Conversely, 
experimental methods involving agglomeration at high temperature must be 
improved to minimize or eliminate effects of backmixing and grid impac- 
tion, and results should be verified by light scattering measurements. 
Experimental results suggest that fractal growth and structure are 
experienced by many agglomeration processes governed by the Brownian 
motion of partrcles*and? clusters? Pit’ is *uniikely that’D,- isa universal 
constant, but that differences in the environment (e.g. T, P, particle 
size and concentration) play an important role in cluster growth and 


SELUCCUTEe: 
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VIII. Recommendations 


The objective of this study was to characterize the structure of 
soot agglomerates in terms of the fractal dimension D, from two dimensio- 
nal projections of three-dimensional shapes. It is important that our 
findings be verified by methods that enable measurements in three 
dimensions. A promising approach is the reconstruction of the 3-D 
structure given projections at several tilt angles. At least two 
reference points, perhaps in the form of surrounding EG ree are 
necessary to insure the proper accounting of particles. Another method 
is to measure the structure factor S(q) for scattering which, from a4 


fractal, is expected to behave as 


S(qyaenqint (11) 


over some range of the momentum transfer q [Weitz and Huang,1984]. 


This study can also be extended to the examination of the reletion- 
ships between the geometric size (R, or LW'/*) and shape L/W to other 
measurable properties such as aerodynamic diameters and diffusion 
constants. In the latter case, the applicability of the Stokes-Einstein 


relation 
D = kT/f (26) 
where D = diffusion constant 


k = Boltzmann’s constant 
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T = absolute temperature 
f = friction force, equal to 3rpd., in the continuum regime, 
where » is the gas viscosity and d, is the particle 


diameter. 


to agglomeration systems can be investigated by expressing d, in terms of 


R, or, in terms of the fractal relation 


d, cea es ye Ps (279 


where N is the number of primary particles in the cluster. 
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Appendices 


Successive Squares Method Program 


PROGRAM NEST 


THIS PROGRAM USES THE SUCCESSIVE SQUARES METHOD WHICH 
OUTPUTS N VS L(LENTH OF SIDE). 

A SERIES OF CONCENTRIC SQUARES OF SIDE L ARE GENERATED 
ABOUT THE CENTER OF MASS AND THE NUMBER OF OCCUPIED 
SITES WITHIN EACH SQUARE IS COUNTED. 


= 


DEFINITION OF VARIABLES 


PIXD(X,Y)= 512X512 MATRIX CONTAINING 1'’S(OCCUPIED) 
AND O’SCUNOCCUPIED). 

Cz, Cx CENTER OF MASS IN TWO DIMENSIONS 

ZX,YY = SUM OF ALL VALUES OF X,Y COORDINATES 

ios eNGIH OF SIDE . 


TNTEGERMPIX DOSS. 51 2)COUNT DE. DY).CX2CY* SIDE “LS, HS 
REAL Xeeyy> >. 

OPEN(1,FILE='’INPUT’ ) 

OPEN(2,FILE=‘OUTPUT’ ) 


et ee OEP RPIXDSARRAY 9 * 45% ** 


DO 10 =rerAs ve 
DOL 20 -J=1 512 
PIZD(T yo )=0 
CONTINUE 
CONTINUE 


***s*** READ DATA AND COMPUTE CENTER OF MASS TO NEAREST INTEGER 


XX=0 

YY=0 

FEAL ea 

DOD SOlL=1), NP 
READ(C) "Dx DY 
PIZD( DX, DY)=1 
XX=XX+DX*1.0 
eet Dew be 


fr 


50 


Gm a 


018) 


() 
O O 


AAAAAL 


CONTINUE 
XX=XX/(NP*1.0) 
YY=YY/(NP*1.0) 

CX=INT(XX+.5) 

CY=INTCYY+ 5) 


* KK KK KK MAIN LOOP kK KK KK K 


COUNT=0 
IF CPIZDC CL CrVEQOal) COUNT=1 
LS=400 
HS=LS/e 
DO 40 I=15Hs 
DOVSO R=— tele 
IF (PIXD(CK+K,CY+I). 
CONTINUE 
DOA GO Rai 
TF ECPEADOCE <1 -CXYrR 
CONTINUE 
DO GOK sate ei 1 eal 
TPeGPI RD CCE+K SCY 1s 
CONTINUE 
DO SO) k=—Ie oli 
TF CPIZD(CZ—-IeCyY+E De 
CONTINUE 
SIDE=1*2+1 
WRITECO Wc DER COUNT 
CONTINUE 


*xexxk*ex PITNTSH NESTED SQUARES 


STOP 
END 
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EQ.1) COUNT=COUNT+1 
EQ.1) COUNT=COUNT+1 
EQ.1) COUNT=COUNT+1 


EQ.1) COUNT=COUNT+1 


PROCEDURES ee. 


Be 


CJC CCl Aca Ae Cr.C2 Co C202 C2 C) Cs CO .02 29-0 C2 09 C2 CD 


C202... 


OQ 


Pair Correlation Method Program 


PROGRAM PAIRCO 


PROGRAM TO COMPUTE C(R) USING MOUNTAIN’S METHOD 
READS IN 2-D SIMULATION OR EXPERIMENTAL DATA ; 
LENGTH BETWEEN ANY TWO PARTICLES ARE COUNTED ONLY ONCE. 


DEFINITION OF VARIABLES 


DISTANCE BETWEEN TWO PARTICLES OR OCCUPIED SITES 


DIST = 
Downe =U LOL gece 
N(I) = NUMBER OF TIMES THE INTEGER DISTANCE I OCCURED 


BETWEEN ALL OCCUPIED SITES. 
PCORR = ACTUAL VALUE OF THE PAIR CORRELATION FUNCTION 
OBTAINED FROM THE FOLLOWING EQUATION : 


PCORR = (2.0*N(1I))/(NP*MAX(I)) 


NP = TOTAL NUMBER OF OCCUPIED SITES 

MAX(I) = TOTAL NUMBER OF POINTS IN A SHELL OF INNER 
RADIUS I-.5 AND OUTER RADIUS I+.5 IN A SQUARE 
LATTICE. 

ARRAY OF X COORDINATES 

ARRAY OF Y COORDINATES 


5d bd 
sa 


REAL X( 156000) ,Y(156000) , DIST, LSQR, PAIRCO 
INTEGER N( 1000), LENGTH 


#4xs**+ OPEN STATEMENTS ******* 


OPEN(UNIT=3,FILE='INPUT'’ ) 
OPEN( UNIT=4,FILE='OUTPUT’ ) 


a 


aqaagaAgd 


CVGy Cease 


20 


aANRA 


40 
10) 


QAaAaAaANg 


5O 


CHG) C25 CG) 


s****** ZERO LENGTH ARRAY ******* 


DO 10 I=1,1000 
N(I)=0 
CONTINUE 


KK RR ROKK 


+xaxx%* READ DATA FROM FILE. 


READ(3,*) NN 

DO 20 I=1,NN | 
READ(S. *) S@lD eee 

CONTINUE 


% eee SOBEGIN NESTED LOOPSTOCCOMPUTE? NCI) S.* 7s 


DO 30 I=1,NN-1 
DO 40 J=I+1,NN 
LSQR=(X(J)-X(1I))**2 +(Y(J)-Y(I))**2 
DIST=LSQR**.5 
LENGTH=INT(DIST+.5) 
N( LENGTH )=N( LENGTH )-+1 
CONTINUE 
CONTINUE 


*4eeee% OUTPUT OF RESULTS ******* 


DO 50 I=1,400 
DRO CNGI ): GIO IeWwRITH ( aete Te Nat) 
CONTINUE 


RRR KK END OUTPUT * KK KK KK 


SLOP 
END 
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SESE CROLL ORO RS Taek eo 


C2 AANA QAAAANQD 


Cent) Gy Cle 


20 


PROGRAM NMALIZ2 


THIS PROGRAM COMPUTES C(R) FROM THE N(I)’S GIVEN 

BY PROGRAM PAIRCO. THE MAXIMUM NUMBER OF OCCUPIED SITES 
FOR EACH VALUE OF R IS STORED IN FILE ‘DATMX2D’ IN 2-D 
UP TO R=250 AND ‘DATMX3D’' IN 3-D UP TO R=30. 


REAL MAX(250) 


INTEGER I,d,DUMMY 
CHARACTER IN*7,OUT*7 


* KKK KK *K OPEN STATEMENTS * Kk KK KK 


PRINT*, ‘ENTER INPUT FILE,OUTPUT FILE (APOSTROPHES NEEDED): ' 


READ, INFOUT 
PRINT* 4° *ENTER NE oITHE GNUMBER “OF -OCCUPIED SITES’ 


: READ* NP 


mK KK OK K OK OPEN STATEMENTS xk RK KK K 


OPEN(1,FILE=IN) 
OPEN(6 , FILE=OUT) 
OPEN(2,FILE='DATMX2D’' ) 


DO WO gla 10 
READ(2,*) J, DUMMY 
MAX(J)=DUMMY*1.0 

CONTINUE 


fx solve! NORMALIZED PALRUO) “Gee es 


DO’ 20) D=1,250 
READ(1,*) J,PCORR 
PCORR= PCORR*2.0*(1I*1.0)/((NP*1.0)*MAX(I)) 
WRITE(6,*) I,PCORR 

CONTINUE 

STOP 

END 
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Ge 


Diffusion Limited Aggregation Program 


OLORS LORS EO Ree lee retake lore te le kee re reel ue ho ky a 


GeO) CoO). C2: C02) 


PROGRAM MK855 


PROGRAM TO SIMULATE DIFFUSION LIMITED AGGREGATION IN 
THREE DIMENSIONS. 


— 


DEFINITION OF VARIABLES 


X(),Y(),2() = 3-D COORDINATES OF PARTICLES BELONGING 
TO THE CLUSTER. 
STP = DISTANCE IN LATTICE SPACINGS TRAVELED BY THE 
PARTICLE PER STEP 
DX,DY,DZ = CURRENT LOCATION OF DIFFUSING PARTICLE 
RMAX = DISTANCE FROM THE ORIGIN TO THE FARTHEST 
PARTICLE IN THE CLUSTER. 
RSTART = RMAX+5, RADIUS ALONG WHICH PARTICLE BEGINS 


TO; DIFFUSE 
KILL = 2.O0*RSTART, DISTANCE BEYOND WHICH PARTICLE IS 
DESTROYED. 
DISTCHK = DISTANCE BETWEEN TWO PARTICLES TO CHECK FOR 
COLLISION. 


PHI,THETA = SPECIFIES SPHERICAL COORDINATES OF DIFFUSING 
~ PARTICLE STARTING POINT ALONG WITH RSTART 
DIRECT = SPECIFIES DIRECTION RANDOMLY (6 POSSIBLE) 


INTEGER X(10000),Y(10000),Z2(10000),'STP, DX , DY "Dz , DIRECT 
REAL RMAX, DUMMY, RSTART,KILL, DIST, DISTCHE 


OPEN STATEMENT FOR OUTPUT 


OPEN(UNIT=2,FILE='OUTPUT’ ) 


RANDOM NUMBER GENERATOR SEED= ISEED 


Z=1005 
CALL RANSET(Z) 
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AaANANM C1 02): C2 


€2: 62.027 7) 


S00 


C2 0.2-C2 C7,@) 


400 


tte) (2 C2 


INITIALIZE SEED COORDINATES 


X(1)=0 

Y(1)=0 

Z(1)=0 

WRETE COU ) KCl )AYCMINZCH) 
RMAX=0 


BEGIN LOOP FOR DIFFUSION LIMITED AGGREGATION 


DO 10 I=2,5000 

DOA=ACT-15."1'.0 

BBY=Y(T-1)7 1.0 

DDZ=Z(I-1)*1.0 

DUMMy, =CDDX *2+DDY" 2+DD24*2)%*2 5 
IF(DUMMY.GT.RMAX) RMAX=DUMMY 
RSTART = RMAX+5.0 


a CHANGED KILL FROM 2*RSTART TO-1.5*RSTARTU******** 


~ 


KILL = 2.0*RSTART 
CONTINUE 


RANDOM NUMBER FOR THETA AND PHI - SPHERICAL COORDINATES 


THETA=RANF()*3.141592654 
PHI=RANF( )*2.0*3.141592654 
DX = INT(RSTART*SIN( THETA) *COS( PHI) ) 


DY = INT(RSTART*SIN( THETA) *SIN( PHI) ) 
DZ = INT(RSTART*COS( THETA) ) 

CONTINUE 

STP = 1 


DIST. = (DES 24DW".2+DZ4 2 Leer 5 


CHANGE STEPSIZE AS THE PARTICLE WALKS IN DIRECTION AWAY 
FROM CLUSTER. 
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IF(DIST.LE.RMAX+10) GOTO 200 
IF(DIST.GT.RMAX+10) STP=2 
IF(DIST.GT.RMAK+20) STP=4 
IF(DIST.GT.RMAX+40) STP=8 
IF(DIST.GT.RMAX+80) STP=16 
200 CONTINUE 
DIRECT=INT(6*RANF( )+1) 


CHOOSE DIRECTION OF WALK RANDOMLY - 6 POSSIBLE DIRECTIONS 


qQaaqaqda 


IF(DIRECT.EQ.1) DX=DX-STP 
IF(DIRECT.EQ.2) DX=DX+STP 
IF(DIRECT.EQ.3) DY=DY-STP 


IF(DIRECT.EQ.4) DY=DY+STP 
IF(DIRECT.EQ.5) DZ=DZ-STP 
IF(DIRECT.EQ.6) DZ=DZ+STP 
DIST DE “825 DY? * 24 D2" 2 tab 
IF(DIST.GT.KILL) GOTO 300 
IF(DIST.GT.RMAX+2) GOTO 400 


DETERMINE WHETHER A HIT OCCURS. 
HIT IF DISTANCE TO ANY PARTICLE IN CLUSTER = 1. 


AQAaAAQARAAaANQ 


DOR204 0-21-12 
DISTCHK=(X(d)-DX)**2 + (Y(d)-DY)**2 + (200 =e 
IF(DISTCHK.EQ@.1) GOTO 40 
20 CONTINUE 


REPEAT “RANDOM WALK EPANO HIT: ELSE, GO 10 40, 


QO'aQAaAa 


GOTO 400 
40 CONTINUE 

X(I)=DX 

Y( I) =DY 

ZCLJ#DzZ 

WRITE(2.79. XC1) YY CE aes) 
10 CONTINUE 


END DLA LOOP. 


AAAANA 


END 


82 


D. Monte-Carlo Digitization Program 


PROGRAM DIGITIZ 


AANAAAAANAAY 


C 
G PROGRAM TO ‘DIGITIZE’ THE CLUSTER GENERATED FROM MOUNTAIN’S 
c PROGRAM TO SIMULATE THE ACTUAL DIGITIZATION PROCESS. RESULTS 
oS WILL BE USED TO EVALUATE THE PAIR CORRELATION FUNCTION. 
#3 
REAL X(9999),Y(9999) ,2( 9999) 
OPEN(UNIT=1,FILE='‘INPUT’ ) 
OPEN(UNIT=2,FILE=‘OUTPUT’ ) 
Z=1931 
CALL RANSET(Z)~ 
CG 
‘C READ COORDINATES FROM FILE 
C 
READ(1,*) NP 
DO 10 I=1,NP 
READ@12 0) SClor mul YL) 
0 CONTINUE 
BEGIN MONTE-CARLO INTEGRATION PROCEDURE 
FOR FILE DATAABQ@, XMIN=11 XMAX= 59 
YMIN=-49 YMAX=18.2 
ZMIN=-4 ZMAX=47 
DO 20 I=-50,20 _ 
DO 30 J=-5,49 
. COUNT=0 
DO 40 K=1,100 
RANX=RANF()-0.5 
RANY=RANF()-0.5 
TX=(I*1.0)+RANX 
TY=(J*1.0)+RANY 
DO 50 L=1,NP 
Dit =) Gite (Lope eal TYR-V CLD) 2% 5 
IF(DIST.LE.0O.5) THEN 
COUNT=COUNT+1 
GOTO 60 
END IF 
50 CONTINUE 
60 CONTINUE 
IF(COUNT.GT.50) THEN 
WRITECG TL vo 
GOTO 70 
END IF 
40 CONTINUE 
70 CONTINUE 


30 CONTINUE 

20 CONTINUE 
STOP 83 
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